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1 Introduction 

The analytic description of wave beam (in particular optical beam) propagation is of 
high importance, both in theory and in practice. 

Because of the mathematical difficulties met when looking for exact solutions, often 
the analytic description must be obtained in an approximate way, and several other times 
one must have recourse to a numerical solution of the (differential or integral) propagation 
equations both in their exact or approximate forms. 

A very common approximation is the paraxial onefTj, rather useful for obtaining 
analytic or numerical solutions. For example, it is by such an approximation that one 
obtains the well-known Fresnel diffraction integral[T], which yields accurate (analytic and 
numerical) results for a large part of the proximal field region, as well as in the transition 
region towards the distant field. 

An important analytic solution forwarded by the Fresnel diffraction integral is the 
gaussian beam, while another one is the Bessel-Gauss beam, found by Gori et al. [2] in 1987. 
The latter, endowed with a transverse proffie in which the Bessel function is modulated 
by a gaussian function, can be regarded as an experimentally realizable version of the 
Bessel beam; indeed, the Bessel beam is a quite noticeable exact, non- diffracting solution 
of the wave equation, but is associated with an infinite power flux (through any plane 
orthogonal to the propagation axis), as it also happens, for instance, with plane waves. 

Notwithstanding the fact that some analytic solutions do exist for the Fresnel diffrac- 
tion integral, they are rare, and normally it is necessary to have recourse to numerical 
simulations. This is particularly true when the mentioned integral is adopted for the 
description of beams generated by finite apertures, that is, of beams truncated in space. 

The past attempts at an analytic description of truncated beams were based on a Fres- 
nel integral: probably the best known of them being the Wen and Breazele method[3], 
using superpositions of gaussian beams (with different waist sizes and positions) in order 
to describe axially symmetrical beams truncated by circular apertures. In that approach, 
those authors had to adopt a computational optimization process to get the superposi- 



tion coefficients, and the beam waists and spot positions of the various gaussian beams; 
actually, the necessity of a computational optimization to find out which beam superpo- 
sition be adequate to describe a certain truncated beam is due to the simple fact that the 
gaussian beams do not constitute an orthogonal basis... 

To cope with the difficulty related with such a computational optimization. Ding and 
Zhang[l] modified the method by choosing since the beginning the beam waist values, and 
then writing down a set of linear equations in terms of the gaussian beam superposition 
coefficients. However, in that new approach the nonhomogeneous terms are given by 
integrals that, once more, cannot in general be easily evaluated in closed form. 

In this paper we are going to show that an analytic description of important truncated 
beams can be obtained by means of Bessel-Gauss beam superpositions, whose coefficients 
are got in a simple and direct way, without any need of numerical optimizations or of 
equation system solutions. 

Indeed, our method is capable of yielding analytic solutions for the 3D evolution of 
Bessel beams, plane waves, gaussian beams, Bessel-Gauss beams, even when truncated 
by finite apertures in the Fresnel regime. 



2 The Fresnel diffraction integral and some solutions 

In this paper, for simplicity's sake, we shall leave understood in all solutions the 
harmonic time-dependence term exp{—iut) . 

In the paraxial approximation, an axially symmetric monochromatic wave field can be 
evaluated, knowing its shape on the z = plane, through the Fresnel diffraction integral 
in cylindrical coordinates: 



\l/(p, z) = exp 



^(kz^'^ 



1^ m/'',0) exp (ik^\ J„ (k^\ p'dp' , (1) 



where k = 27r/A is the wavenumber, and A the wavelength. In this equation, p' recalls 
us that the integration is being performed on the plane z = 0; thus, \E'(p', 0) does simply 



indicates the field value on z = 0. 

Let us consider a gaussian behaviour on 2; = 0, that is to say, let us choose the 
"exitation" 



^(p',0) = ^g(p',0) = Aexp{-qp' 



/2^ 



(2) 



contained in ref. [1]; with A and q constants that can have a complex value. For Re(g) > 0, 
one gets the pretty known gaussian beam solution: 
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where 



Q = q — ik/2z 



(4) 



Another important solution is obtained by considering on the z = plane the excita- 
tion given by 



vi>(p',0) = ^bg{p,0) = AJo{kpp')exp{-qp") , 
which, according to eq.Q, produces the so-called Bessel-Gauss beam[2]: 
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quantity Q being given by eq.Q, and kp being a constantF] 

The Bessel-Gauss beam given by eq.(|6| is particularly interesting since it can well be 
regarded as a realistic version (experimentally speaking) of the ideal Bessel beam: 



^b(p,2;) = AJq (kpp) exp{ik^z) 



(7) 



where kz = Jtu'^/c'^ — k'l. 



* Quantity kp is the transverse wavenumber associated with a Bessel beam transversahy modulated by 
the gaussian function. 



The Bessel beam, eq.Q, is an exact solution to the wave equation and is known to 
possess the important characteristic of keeping its transverse behavior unchanged while 
propagating, so to belong to the class of the nondiffracting beams. However, the ideal 
Bessel beam is endowed with an infinite power flux, and cannot be concretely generated. 
By contrast, the Bessel-Gauss beam, eq.([6]), modulates in space the transverse behaviour 
of the Bessel beam by a gaussian function, getting a finite power flux. The Bessel-Gauss 
beam will no longer remain indefinitely undistorted, but nevertheless shows to possess a 
rather good resistance to diffract ion [2]. 

The gaussian beam, eq.([3]), and the Bessel-Gauss, eq.(|6]), solutions are among the few 
solutions to the Fresnel diffraction integral that can be got analytically. The situation 
gets much more complicated, however, when facing beams truncated in space by finite 
circular apertures: For instance, a gaussian beam, or a Bessel beam, or a Bessel-Gauss 
beam, truncated via an aperture with radius R. In this case, the upper limit of the 
integral in eq. ([I]) becomes the aperture radius, and the analytic integration becomes very 
difficult, requiring recourse to lengthy numerical calculations. 

As we already mentioned, in their alternative attempt at describing truncated beams. 
Wen and Breazele|3] adopted superpositions of gaussian beams. To be, now, more specific, 
those authors wrote down the solution for a wave equation in the paraxial approximation 
as: 
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with 
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Solution (|8j) is a superposition of A^ gaussian beams. Its coefficients An and g„ are to 
be obtained starting from the field existing on the z = plane, that is, starting from the 
initial excitation that we shall call V{p). One therefore looks for \l'(p, 0) = V{p). From 
eq.(|8|) one gets 



N 



V{p) = ^A„exp(-g„p2) . (10) 

n=l 

The initial field V{p) can represent a beam (e.g., a plane wave, a gaussian beam,...) 
truncated by a circular aperture with radius R. To get the coefficients An and g„ from 



eq.(lO), those authors had recourse to a computational optimization process in order to 



minimize their mean square error: Namely, to minimize the difference between the desired 



function V{p) and the gaussian series in the r.h.s. of eq.(lO). Such a method yields good 
results, provided that the exitation function V{p) does not oscillate too much. But the 
coefficients A^ and g„ are obtained in a strictly not algebraic manner, depending on the 
contrary on numerical calculations. 

Let us be more specific also about the modification of Wen and Breazele's method 
introduced by Ding and Zhang[4J. They postulated the values of the parameters g^ and, 
by minimizing the mean square error between the desired function V{p) and the gaussian 



series in eq.(lO), arrived at a system of linear equations containing the unknowns An (the 
coefficients of the gaussian beam superposition), without needing a numerical optimization 
process. However, in the system of equations needed to determine the coefficients An, 
the non-homogeneous terms consist in integrals that, depending on the field one wishes 
to truncate (i.e., depending on V{p)) can be difficult to be calculated analytically... 

In the next Section we are going to propose a method, for the description of truncated 
beams, that appears to be noticeable for his simplicity and, in most cases, for its total 
analiticity. Our method is based on Bessel-Gauss beam superpositions, whose coefficients 
can be directly evaluated without any need of computational optimizations or of solving 
any coupled equation systems. 



3 The method 

Let us start with the Bessel-Gauss beam solution, eq.([6]), and consider the solution 
given by the following superposition of such beams: 
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quantities An being constants, and Qn being given by eq.(|9]), so that Qn = Qn — ik/2z, 
where g„ are constants that can have complex values. Notice that in this superposition 
all beams possess the same value of kp. 

Let us recall, incidentally, that all the beams we are considering in this work are 
important particular cases of the so-called Localized Waves (see [HI ?,[7] and refs. therein; 
see also ^^). 



Our purpose is that solution (11) be able to represent beams truncated by circular 



apertures: As announced, we are particularly interested in the analytic description of 
truncated beams of Bessel, Bessel-Gauss, gaussian and plane wave types. 

Given one of such beams truncated at^; = by an aperture with radius R, we have 



to determine the coefficients An and g„ in such a way that eq.(ll) represents with fidelity 



the resulting beam. If the truncated beam on the z = plane is given by V{p), we have 
to obtain \l/(p, 0) = V{p)] that is to say: 
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The r.h.s. of this equation is nothing but a superposition of Bessel-Gauss beams, all with 
the same value kp, at z = [namely, each one of such beams is written at ^ = according 
toeqM]. 



Equation (12) will provide us with the values of the A^ and g„, as well as of A^. Once 



these values have been obtained, the field emanated by the finite circular aperture located 



at 2; = will be given by eq.(ll). Remembering that the g„ can be complex, let us make 
the following choices: 



Qn = QR + iqir. 



(13) 



where q^ > Q is the real part of g„, having the same value for every n, and g/.„ is the 



imaginary part of g„ given by: 



Qln = -^^' (14) 



where L is a constant with the dimensions of a square length. 



With such choices, and assuming A^ — )■ oo, equation (12) gets written as 

^(p) = MkpP) axp (-qRp'^) J2 Anexp(i——p^] , (15) 

n=-oo \ L / 

which has then to be exploited for obtaining the values of An, kp, q^ and L. 

Let us recall that aim of our method is describing some important truncated beams 
starting from the value of their near fields (i.e., of their fields in the Fresnel region). 

In the cases of a truncated Bessel beam (TB) or of a truncated Bessel-Gauss beam 



(TBG), it results natural to choose quantity kp in eq.(15) to be equal to the corresponding 
beam transverse wavenumber. 

In the case of a truncated gaussian beam (TG) or of a truncated plane wave (TP), 



by contrast, it is natural to choose kp = in eq.(15). 
In all cases, the product 



exp(-g/?p^) Yl ^nexpU-— p^j , (16) 



in eq.(15) must represent: 

(i) a function circ(p/-R), in the TB or TP cases; 

(ii) a function exp(— gp^) circ(p/i?), that is, a circ function multiplied by a gaussian 

function, in the TBG or TG cases. Of course (i) is a particular case of (ii) with g = 0. 

It may be useful to recall that the circ-function is the step-function in the cylindrically 

symmetrical case. Quantity R is still the aperture radius, and circ(p/i?) = 1 when 

< p < R, and equals in the contrary case. 



Let us now show how expression ( 16 ) can approximately represent the above functions, 
given in (i) and (ii). 



To such an aim, let us consider a function G{r) defined on an interval \r\ < L/2 and 
possessing the Fourier expansion: 



G{r) = Yl Anexp{i27mr/L) for \r\ < L/2 



(17) 



where r and L, having the dimensions of a square length, will be expressed in square 
meters (m^). 

Suppose now the function G{r) to be given by 



G{r) 



exp (g/jr) exp (— gr) for \r\ < R"^ 
for i?2 < \r\ < L/2 , 



where g is a given constant. 

In this case, the coefficients An in the Fourier expansion of de G{r) will be given by: 
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in eqs.(T7p8) we shall have that 
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exp(gi:jp^) exp(— gp^) for |p| <R 



for i? < IpI < J L/2 



(21) 



where the coefficients An are still given by eq.(19). One recognizes that the l.h.s. of 



eq.(21) is the term multiplying the gaussian in expressio (16). 

The l.h.s. of eq.(21), which depends on p^, will not be exactly periodical: But it will 



re-assume the values, assumed in the fundamental interval (|p| < w(L/2)), in shorter and 



shorter further intervals, with decreasing spatial "periodicity", for \p\ > J(L/2). In such 



a way, with p > 0, expression (16) becomes: 
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where /(p) is a function existing on decreasing space intervals, and assuming exp[(g/j — 



q)R'^] as its maximum values. Since \/L/2 > R, for suitable choices of qn e L, we shall 



have that exp(— g^jp^) /(p) ^0 for p > JL/2. 
Therefore, we get that 



expi^-qup'^j Y^ Anexp{i27mp'^/L) ^ exp(^-gp^j circ(p//?) , (23) 
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which corresponds to the case (i), when g = 0, and to the case (ii). 



Let us recall once more that the An are given by eqs.(19). 



On the basis of what shown before, we have now in our hands a very efficient method 
for describing important beams, truncated by finite apertures: Namely, the TB, TG, 
TBG, and TP beams. Indeed, it is enough to choose the desired field, truncated by a 



circular aperture with radius R, and describe it at z = by our eq.(15). Precisely: 



In the TBG case: the value of kp in eq.(15) is the transverse wavenumber of the 



Bessel beam modulated by the gaussian function; A„ is given in eq. ( 19 ) ; g is related 



to the gaussian function width at z = 0. The value L and qn, and the number N of 



terms in the series (15), are chosen so as to guarantee a faithful description of the 



beam at 2; = when truncated by a circular aperture with radius R. 

In the TB case: the procedure will be the same as for TBG, but with q = 

In the TG case: the procedure is still the same as for TBG, but with kp = 



10 



• In the TP case: the procedure is once more similar, but using this time kp = e 
q = 

Finally, once the chosen beam is described on the truncation plane (z = 0), the beam 



emanated by the finite aperture will be given by solution (11). 

It is important to notice at this point that, for a given truncation, innumerable sets 
of values of qr and L exist which yield a faithful description of the truncated field. The 



choice is made in order that the solution given by the series (15) has good convergence 
properties. Of course, one will use a finite number 2N + 1 of terms in the mentioned 
series, with —N < n < N. 

Let us go on to some examples. 

4 Applying the method 

In this Section we shall apply our method, as we already said, to situations in which 
important truncated beams appear. Notice that below we shall always assume a wave- 
length of 632.8 nm. 

4.1 Analytic description of the truncated Bessel beam 

Let us start from a Bessel beam, truncated at 2; = by a circular aperture with radius 
R; that is to say, from \E'7^b(p, 0) = Jo{kpp) ciYc{p/R). 

Let us choose R = 3.5 mm, and the transverse wavenumber kp = 4.07 ■ 10^ m~^, 
which corresponds to a beam spot with radius approximatively equal to Ap = 59 fim 
(while A = 632.8 nm, as always). 



At 2; = the field is described by eq. ( 15 ) , where the An are given by eq. ( 19 ) and where 
g = 0. In this case, a quite good result can be got by the choice L = 3R^, qn = Q/L and 
N = 23. Let us repeat that, since such a choice is not unique, very many alternative set 
of values L and qr exist, which yield as well excellent results. 



ure 1 shows the field given by eq.(15): it represents with high fidelity the Bessel beam 
truncated at z = 0. 
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Figure 1: Field given by eq.(15), representing a Bessel beam at z = 0, with kp = 4.07 



10^ m ^ and truncated by a finite circular aperture with radius R = 3.5 mm. The 
coefficients An are given by eq.(19), with q = 0, L = 3R^, qr = 6/L and N = 23. 



The resulting field, emanated by the aperture, is given by solution (11), and its 



intensity is shown in Fig. 2. One can see that the result really corresponds to a Bessel 
beam truncated by a finite aperture. Figure 3 depicts the orthogonal projection of the 
same result. 



Increasing A^, that is, increasing the number of terms in the series ( 11 ) which expresses 



the resulting field, while keeping the same values for L and qn, the spatial shape of the 
obtained field practically will not change; but there will become more evident the rapid 
oscillations that occur at the beam crest, i.e., for '^tb{p = 0,-2). This is shown in Fig. 4, 
were we used N = 500. 

4.2 Analytic description of the truncated gaussian beam 

Let us go on now to consider a gaussian beam truncated at 2; = 0; that is, \E'tg(Pj 0) = 
exp(— gp^) circ(p/i?), whose initial intensity spot radius is Ap = 59 /im, and therefore 
q = l/(2Ap^) = 144 • 10^ m^^. The radius of the circular aperture is equal to the beam 
spot radius, i.e., i? = 59 /im, while, as always, A = 632.8 nm. 



The situation at 2; = is still described by eq.(15), with kp = 0, where the An are 



given by eq.(19). Now, a good result can be obtained, for instance, by using the values 
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Figure 2: Intensity of a Bessel beam truncated by a finite aperture, as given by solution 



(11). 




Figure 3: Orthogonal projection of the intensity shown in Fig. 2. 



L = iR^, qr = 8/L and A^ = 81. 



In Fig.5 we show the field given by eq.(15), in the case of a gaussian beam truncated 



at z = 0. The dotted line depicts the ideal gaussian curve, without truncation. 



The resulting field, emanated by the finite aperture, is given by the solution (11) 



and Fig. 6 shows its square magnitude. In Fig. 7 we show the corresponding orthogonal 
projection. 
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Figure 4: Oscillations of the field intensity on the z axis, when adopting A^ = 500 in 



eq.(ll). 
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Figure 5: In this figure it is shown the field given by eq.(15), when representing a gaussian 
beam truncated at 2; = 0. The dotted line depicts the ideal gaussian curve, when 
truncation is absent. 



4.3 Analytic description of a truncated Bessel-Gauss beam 



Let us now consider the interesting case of a Bessel-Gauss beam truncated by a 
circular aperture of radius R\ that is, \I/t_bg(p?0) = jQ{kpp)ey.'p{—qp'^)c\ic{p/R), with 



fcp = 4.07-10^m-\ q 



1.44 ■ lO'^m" 



R = 1 mm, and A = 632.8 nm. 
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Figure 6: This Figure shows the square magnitude of the field emanated by a finite 



aperture, in the case of a gaussian beam, according to the solution (11). 
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Figure 7: This Figure shows the orthogonal projection coresponding to the previous 
Figure. 



The situation at z = is described by eq.(15), where the An are given by relations 
(19). A very good result can be otained, e.g., by adopting the values L = 10 R^, Qr = q 

15 



and A^ = 30. Figure 8 shows the field in eq.(15), in the present case of a Bessel-Gauss 
beam truncated at z = 0. 



The resulting field emanated by the finite aperture is given by solution ( 11 ), and Fig. 9 



shows its square magnitude. Figure 10 depicts the orthogonal projection for this case. 
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Figure 8: This Figure shows the field in eq.(15), in the case now of a Bessel-Gauss beam 
truncated at ^ = 0. Here we adopted the values L = 10 R^, Qr = q and A^ = 30. See the 
text. 
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Figure 9: This Figure shows the square magnitude of the field emanated by a finite 



aperture in the case of a truncated Bessel-Gauss beam, represented by solution (11) 
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Figure 10: This Figure depicts the orthogonal projection for the case in the previous 
Figure. 
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4.4 Analytic solution of a truncated plane wave 

Consider now the case of a plane wave truncated by a circular aperture em z = 0, 
that is, \l/j'p(p, 0) = circ(p/i?), were we choose R = 1 mm and A = 632.8 nm. 



Once more, eq.(15) describes the field at z = 0, with kp = 0, the coefficients An being 



given by relations (19), with g = 0. A good result can be got adopting, e.g., the values 



L = 6R^, qR = 8/L and A^ = 150. 



Figure 11 shows the field in eq.(15), in the present case of a plane wave truncated at 
z = 0. 



The resulting field emanated by the finite aperture is given by solution (11), and its 



square magnitude is shown in Fig. 12. Figure 13 depicts the corresponding orthogonal 
projection. 
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Figure 11: This Figure shows the field at z = 0, with kp = 0, as given by eq.(15 in the 
new case of a plane wave truncated at 2; = 0. We have here adopted the values L = 6R^, 



Qr = 8/L and A^ = 150. the coefficients An being given by relations (19), with g = 0. 
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Figure 12: This Figure shows the square magnitude of the resulting field, emanated by 



the finite aperture, as given by solution ( 11 ) in the present case of a truncated plane wave. 
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Figure 13: The orthogonal projection corresponding to the previous Figure. 

5 Conclusions 

In this paper, starting from suitable superpositions of Bessel-Gauss beams[2], we have 
constructed a simple, effective method for the analytic description, in the Fresnel region. 
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of important beams truncated by finite apertures. 

The solutions obtained by our method, and representing truncated Bessel beams, trun- 
cated gaussian beams, truncated Bessel-Gauss beams and truncated plane waves, fully 
agree with the known results obtained by lengthy numerical evaluations of the correspond- 
ing Fresnel diffraction integrals. (Incidentally, let us mention that all the beams considered 
in this work are important particular cases of the so-called Localized Waves [S, ?, [TJ E])- 

At variance with the previous Wen and Breazele's approach|3] (which uses a computa- 
tional method of numerical optimization to obtain gaussian beam superpositions describ- 
ing truncated beams), and even at variance with Ding and Zhang's approach[3j (which is 
an improved version of Ref.[3]), our method does not need any numerical optimizations, 
nor the numerical solution of any coupled equation systems. 

Indeed, the simpler method exploited in this paper is totally analytic, and directly 
applies to the beams considered above, as well as to many other beams that are being in- 
vestigated and will be presented elsewhere: like truncated higher order Bessel and Bessel- 
Gauss beams; or beams truncated by circular apertures; or beams truncated and modu- 
lated by convergent /divergent lenses; etc. In particular we have applied this method to 
remote sensing by microwaves (cf.,e.g., Ref. [10]), constructing finite antennas which emit 
truncated Bessel beams with the required characteristics (patent pending). Of course, 
this method works in Electromagnetism (Optics, Microwaves,...), as well as in Acoustics. 

Let us stress that one of the main byproducts of our mathematical method is that 
by it one can get in few seconds, or minutes, high-precision results which could otherwise 
require several hours, or days, of numerical simulation. 
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